Fast migration and emergent population dynamics 
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We consider population dynamics on a network of patches, each of which has a the same local 
dynamics, with different population scales (carrying capacities). It is reasonable to assume that if the 
patches are coupled by very fast migration the whole system will look like an individual patch with 
a large effective carrying capacity. This is called a "well-mixed" system. We show that, in general, 
it is not true that the well-mixed system has the same dynamics as each local patch. Different global 
dynamics can emerge from coupling, and usually must be figured out for each individual case. We 
give a general condition which must be satisfied for well-mixed systems to have the same dynamics 
as the constituent patches. 
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Many population models consider patches of habi- 
tat, each with a local birth-death dynamics, where the 
patches are coupled by diffusive migration [l|, Q . Here we 
will consider a special case where all the patches have the 
same local dynamics except for a parameter that gives 
the scale of the local mean population (carrying capac- 
ity). That is, we consider a metapopulation which can 
exchange between patches which are more or less favor- 
able places to be. 

Such models are notoriously complex and have rich 
behavior However, one particular case seems simple: 
suppose the migration rate is very fast, faster than any 
other rate in the problem. In this case all the populations 
will become equal on average. It is tempting to think that 
the autonomous dynamics of the total population will be 
the same as that on each patch with some sort of average 
carrying capacity, i.e., the whole system will act as one 
patch carrying the total population [3j|. In this case the 
system is "well-mixed" . 

In this paper we will show that this common belief 
is false in general: there are tight restrictions for the 
type of dynamics that admits this type of homogeniza- 
tion. If these restrictions are not obeyed the dynamics 
of the emergent behavior of the total population can be 
different in nature from that of the individual patches. 
We will derive the restrictions on the dynamics and give 
examples of the sort of emergent behavior that we have 
described, first by using mean-field theory and then in- 
cluding fluctuations. 

We first formulate in a more precise form the meaning 
of the idea that a total population of a number of local 
patches evolves like a larger population on a single patch. 
Let X be the size of a population. The population dy- 
namics is driven by the birth-death processes X — > X+r, 
each process associated with the rate W(X;r). We shall 
assume that the rate can be put into the form 

W(X;r) = Kw(x;r) 
x = K~ X X. 

The number K is the carrying capacity. Now consider 



a connected network of N patches such that the local 
dynamics of the population Xi is identical on each patch, 
but the carrying capacities, Ki are different. That means 
the local birth-death processes are of the form Xi — > 
Xi + r, with corresponding rates 



Wi(Xi; r) = Kiw(xi;r), i 
x t = K- l X t . 



1. 



(1) 



The patches are connected by migration. We define 
the dynamics of the total population X — J2i x i to be 
well-mixed if it is driven by transitions X — > X + r, with 
corresponding rates 



W{X; 



Kw(x;r) x = K~ 1 X, 



where w[x;r) are local rates (|T|). The number K is the 
effective carrying capacity of the total population. This 
defines well-mixed population dynamics. 

Now assume that the rate of migration between 
patches is much larger than the local birth-death rates. 
We first consider the mean field limit, neglecting fluc- 
tuations which is appropriate for Ki ~ K ^> 1. The 
mean-field equations are 

X t = J2 r Wi(Xi;r) + Dj2( x j-Xi) 

r jeli 

= ]T rKiwiK^Xf, r) + D^2 i X o ~ X i) ( 2 ) 

where D is migration rate, which is taken to be identical 
for all groups, and Xi is the set of indices, associated with 
the patches, connected by migration to patch i. We are 
interested to see if for D — > oo, the total population is 
driven by the homogenized dynamics: 



X = ^rKw{K~ 1 X;r), 



(3) 



for some effective carrying capacity K. 

Consider a specific example first where the local dy- 
namics displays the Allee effect i.e. in a certain range 



2 



of parameters the single patch dynamics is bistable: there 
is one stable state with finite population and another cor- 
responding to local extinction. Our example for the local 
birth-death process is: 

X t — > 0; Rate [iX t 
2Xi -» 'iXi- Rate (\/2K i )X i (X l - 1) 
3Xi -»• 2X,- Rate ( f j/6Kf)X l (X l - l)(X l - 2), 

on each site, where Ki is a local carrying capacity. We 
note that the rates have the canonical form ((T|). The 
mean-field equations are: 



X t 



A 2 a 



6Kf 



■xt 



(4) 



For D = the local stable stationary states are found to 
beX* = {0,^f i [l + (l-|^) 1 / 2 ]}, i.e., local bistability 
(Allee effect) is observed for A 2 > (8/3)ct/x. 

For large migration rates the local populations are ap- 
proximately equal, Xi = X/N, on the time-scale of lo- 
cal dynamics. The mean-field equation of motion for 
the total population can be derived by summing up the 
equations Q to get rid of migration terms, and setting 
X, = X/N: 



X 



+ Xx 2 - ^Lx 3 

2K 6K 2 
/ \ -i 



(5) 



Comparing Eqs.([5]) and ((3]) we see that the total popula- 
tion X follows a single-patch dynamics with the rescaled 
rate a. Note that for equal local carrying capacities, 
Ki = Kj for all a = a. On the other hand, for 
different carrying capacities it can be shown that er > a, 
i.e., the well-mixing condition is not satisfied. 

This can lead to a qualitative difference of the local 
dynamics and dynamics of the total population: specifi- 
cally, the local dynamics is bistable for A 2 > |cr/i, while 
for A 2 < |ajU the only stable state is extinction of the 
population. Since a > a, in a certain range of parameters 
one can simultaneously have A 2 > |<7/Lt and A 2 < In 
this case, even though each local population separately 
can be bistable, the total population in the presence of 
fast migration is driven to extinction for any initial distri- 
bution of local populations. Next we present an example 
illustrating this remarkable dynamics. 

We consider the following dependence of the local car- 
rying capacity on the (discrete) spatial coordinate i: 

Ki=K(l + esm(2m/L)), 

so that the carrying capacity is a periodic function of i, 
and the parameter e determines the relative magnitude 



of these spatial variations, see Fig.QJi (dashed lines). For 
sufficiently small e the nontrivial nonzero state is sta- 
ble. However, when e exceeds the threshold determined 
by A 2 = |cr(e)/i, the system is rapidly driven to extinc- 
tion. Fig. QJ, shows the spatial population profiles Xi for 
two values of e (above and below the threshold) ; Fig. [TJd 
shows the corresponding time dependence of the total 
population X. 
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FIG. 1: (a) Spatial population profiles Xi at intermediate 
time t — 100 for two values of normalized variation of carrying 
capacity: e = 0.45 (upper solid line) and e = 0.55 (lower solid 
line). The dashed lines show the corresponding spatial profiles 
for zero diffusion, (b) The total population X as a function 
of time for e = 0.45 (upper solid line) and e — 0.55 (lower 
solid line). The parameters are: K — 60, D — 300, /i = 0.09, 
A = 0.2, a = 0.15. 

In the general case the necessary condition for the well- 
mixing is obtained from Eq.([2]) and ([3]) following similar 
manipulations. The condition for well-mixing is: 

rK lW {{NKi)- l X; r) = ^ rKw (k^X; r) , 
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for some K. 

We can find the necessary and sufficient condition for 
well mixing beyond mean-field theory. Stochastic popu- 
lation dynamics is driven by the master equation for the 
probability distribution: 

N 

P(X) = ^^W 4 (X l -r;r)P(X-re 1 ) 

i— 1 r 
N 

- ]T]T^(X l ;r) J P(X) + .DWoP(X) (6) 



where Wi(Xi] r) are local transition rates, = 5ij and 
DWo (X, r) is the migration operator: 

1 N 

Wb(X,r)P(X) = -]T^pr j + l)P(X-e l +e J ) 

i=i jeii 

1 N 

+ 2EE( I ' + 1 ) p ( x+e '-^ 

i=l jGXi 
i=l jGXi 

We use on-site dynamics identical on all patches with 
different local carrying capacities Ki We also assume the 
carrying capacities are large, Ki = Kn%, K ~3> 1. We 
seek the dynamics of total population X = J2i -Xj in the 
limit D — y oo. 

For K 3> 1 we can use the eikonal approximation for 
the master equation [Bj-TToj , looking for the probability 
distribution in the form 



P(X,t) = e~ KS ^ 



(7) 



For the quasistationary distribution we obtain Hamilton- 
Jacobi equation for the action S (x) ll| : 



(8) 



P(x,p) = 

p t = diS (x) , 

associated with the following Hamiltonian: 
P(x,p) = PP (x,p) + Pi(x,p) 

N 

Pu(x,p) = J2Y, K i W ( K ^ X i> r )( ePir -V 

i—1 r 

Po(x, P ) = i^^^( e ^-l) 

The Hamiltonian © generates the fluctuational dynam- 
ics of the system, in the sense that the most probable 



path of the system, leading from Xj at t = tj to x/ at tt 
can be shown Q to correspond to Hamiltonian dynamics 
generated by (JSJl- To derive the fluctuational equations 
of motion for the (rescaled by K) total population, we 
do a canonical transformation from x and p to Q and 
Pi,P Q ,i = 2,3,...,N: 



q, 



Pi = Pq Pi = Pc 



Q 



P. 



In the new variables, the Hamiltonian is in the form 
H = DH ({ qi },Q,{p gi }) +H 1 ({q i },Q,{p qi },P Q ) with 
i = 2,3,..., N. 

The total population Q and the conjugate momentum 
Pq are slow variables, evolving on the time-scale unity. 
The variables qi,Pi are fast, evolving on the time-scale 
of migration e = D^ 1 . In zero order in e the evolution 
of the total population is driven by the Hamiltonian: 

H slow (Q, Pq) = Pi ({»(<?)} , Q, 0, Pq) , (9) 
9i(Q) = Q/N, i = 2,3,...,N, 

which reflects the adiabatic slaving of the fast local pop- 
ulations qi , by the slow evolution of the total population 
Q, such that instantaneous sizes of the local populations 
are all equal, due to the fast migration. We outline the 
derivation below. This derivation generalizes the fast mi- 
gration result presented in Ref. \2\ . 

The Hamiltonian equations of motion in new variables 
and after rescaling of time t — > Dt are: 

qi = c)p i Po(q,g,P)+ec)p ! P 1 (q,Q,P,PQ) 

P = -^P (q,Q,P)-ea :c ,P 1 (q,Q,P,PQ) (10) 

Q = £<9pPi(q, Q, P, Pq) 

Pq = -a Q Po(q,Q,P)-ea Q Pi(q,Q,P,PQ), 

where we use notation q = (<?2, 93, <Ziv) and P — 
(P 2 ,P 3 ,...,P N ) for brevity. 

We look for = qi(Q) + eq^' and Pi = eP^, where 
qi{Q) is the stationary value of qi for the instantaneous 
value of Q. In the leading order in e the first two lines in 
equations of motion (| 10[) become 

= d Pl P (q,Q,P) 

= -^Po(q,Q,P)-^ i P 1 (q,Q,P,P Q ) (11) 

The first of Eqs. pTj) determines qi(Q), which is Q/N for 
the Hamiltonian ([5]). In view of the functional form of 
the Hamiltonian, H (q, Q, P) = P • h(q, Q) + 0(e 2 ), we 
obtain from the first of Eqs. (ITTT) 

dH (q,Q,P) dq 3 gg (q,Q,P) _ 

d qj dQ dQ { ' 

to first order in e. Combining the second of Eqs. ([TTj) and 
Eq. ljTlZ)) we derive 

e ggi(q, Q, P, Pq) <% _ dH (q, Q, P) 



dqj 



dQ 



dQ 
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Using Eqs. flT3|) and coming back to the original dimen- 
sional time we rewrite the last two Eqs. (|10p in the form 



Q = d P Hi(q(Q),Q,0,P Q ) = d P H slow (Q,P) 
P = -d x . Hl (<i(Q),Q,0,P Q )^L 

- dQHifriQ^Q^PQ) = -d Q H slow {Q,P), 

where H slow (Q,P Q ) = fli(q(Q), Q, 0, P Q ), Eq.©. 

For simplicity of notation we omit the subscript of Pq 
in what follows. Using Eqs.© and (J9j> we derive: 



JY 



H s i ow (q, p) = y^y^ 



- 1)(14) 



This Hamiltonian drives the fiuctuational dynamics of 
the total population Q. 

In order that the total population follows an effective 
single-site dynamics with a rescaled carrying capacity 
(well-mixing), the Hamiltonian (1141) must have the form 
of a single-site Hamiltonian: 



h (q, p) = J2 kw r ) ( ePr - !) 



(15) 



The equivalence of expressions (|T4"f and (|T5j) gives the 
necessary and sufficient condition for well-mixing: 



A' 



= kw (h 1 Q;r), Vr, (16) 



where k = KK~ X is the rescaled effective carrying capac- 
ity. To make further progress we assume that the on-site 
birth-death rates w {x; r) can be Taylor-expanded around 
x = 0: 



■ 0; r) = < 



Vr, x. 



Then condition (|16|) is equivalent to 



JV 



»=i 



0, Vr, n. 



(17) 



The general solution of this equation for varying carrying 
capacities, Ki ^ Kj for some is 

a n = a lSn,l + a^«<5n,n* Vr 

l 

TV 
i=l 

where n* is a fixed natural number. In other words, a 
rescaled carrying capacity can be found if and only if the 
on-site birth-death rates w (x; r) have the form 



w (x; r) = a\x + a r n »x n Vr. 



(18) 



We note that for identical carrying capacities, Ki = 1 
for all i, the condition (fl~7|) is satisfied for arbitrary a r n 
and K = N, i.e., for identical carrying capacity the total 
population evolves like a single-site population with a 
carrying capacity rescaled by the number of sites N. 

Many population dynamics models have the form 
for n* = 0, 2. For n* — 0, as in Poissonian process [lc 



k = n(n- 1 J2k^J , 



i.e., the carrying capacity is N times the arithmetic mean 
of the local carrying capacities. For n* = 2, as in logistic 
growth [14j, 



i.e., the carrying capacity is N times the harmonic mean 
of the local carrying capacities. 

In this paper we have shown that a notion that seems 
to be completely obvious is false. It is not true, in general, 
that if we mix very fast, a metapopulation with identical 
dynamics on patches of habitat (with different carrying 
capacities) will synchronize to act like a single population 
with the same dynamics. 

The notion is true in very simple special cases: in ef- 
fect, if the local dynamics is characterized by just one 
parameter (the carrying capacity) then fast migration 
can only average things. Most examples that come to 
mind are exactly of this sort. But if, as in the case of the 
Allee effect, there is another relevant parameter, then it 
can change too, and the qualitative behavior of the syn- 
chronized population can be different from that of any 
individual patch. 

For the general case, only if Eq. (IT6l) is obeyed, do we 
have the same dynamics in the well-mixed case as for the 
individual patches. If the rates are polynomials in the 
populations, only if we have linear plus one other power, 
do we conserve the dynamics. 

We thank B. Meerson for many discussions and interest 
in this work. M.K. thanks M. Dykman for discussions 
on the adiabatic approximation. M.K. is grateful to the 
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